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Abstract. In this paper we consider the population of eccentric bina- 
ries with a neutron star and a white dwarf that has been revealed in our 
galaxy in recent years through binary pulsar observations. We apply our 
statistical analysis method (Kim, Kalogera, & Lorimer 2003) and calcu- 
late the Galactic formation rate of these binaries empirically. We then 
compare our results with rate predictions based on binary population syn- 
thesis from various research groups and for various ranges of model input 
parameters. For our reference moel, we find the Galactic formation rate 
of these eccentric systems to be ~ 7 Myr -1 , about an order of magnitude 
smaller than results from binary evolution estimations. However, the em- 
pirical estimates are calculated with no correction for pulsar beaming, 
and therefore they should be taken as lower limits. Despite uncertain- 
ties that exceed an order of magnitude, there is significant overlap of the 
various rate calculations. This consistency lends confidence that our cur- 
rent understanding of the formation of these eccentric NS— WD binaries 
is reasonable. 



1. Introduction 

Binary pulsars with white dwarf companions (NS— WD binaries) in eccentric 
orbits have been revealed with binary pulsar and optical observations in recent 
years. This sub-population is considered to be rather special because, based on 
our standard understanding of binary evolution, NS— WD binaries are expected 
to be circular: a neutron star forms first from the original binary primary and 
the white dwarf formation follows mass transfer episodes that are expected to 
circularize the binary orbit. For a review of this scenario, see Lorimer (2001). 
The neutron star produced by this evolutionary path is expected to be 'recycled', 
i.e. spun-up by mass accretion from its companion. However, the existence of 
eccentric NS— WD binaries such as PSR B2303+46 (Stokes, Taylor, & Dewey 
1985; van Kerkwijk & Kulkarni 1999) and PSR J1141-6545 (Kaspi et al. 2000; 
Bailes et al. 2003) implies that a different evolutionary path from the standard 
scenario is also be possible (Tutukov & Yungelson 1993; Portgies Zwart & Yun- 
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gelson 1999; Tauris k Sennels 2000; Nelemans, Yungelson, k Portegies Zwart 
2001; Brown et al. 2001; Davies, Ritter, k King 2002). The non-zero eccen- 
tricity of the binary orbit is introduced by a supernova (SN) explosion of the 
secondary companion that occurs after the original primary has already evolved 
into a white dwarf. 

Our motivation for this study is to estimate the Galactic formation rate 
of eccentric binaries based on the observed pulsars. We apply our statistical 
analysis developed to estimate the Galactic double-neutron-star (DNS) merger 
rate (Kim, Kalogera, k Lorimer 2003) and derive a probability density func- 
tion (PDF) of the Galactic formation rate for eccentric NS— WD binaries. This 
method provides us with rate estimates that are independent from those ob- 
tained using binary evolution calculations. We compare our empirical rate es- 
timates to results from population synthesis calculations and conclude that, 
despite the large uncertainties, the results are indeed consistent. 



2. Formation of eccentric NS— WD binaries 

In our classical understanding of binary evolution, we expect that the formation 
process of NS— WD binaries in close orbits involves the circularization of the 
binary orbit, even for systems with massive white dwarfs: the neutron star 
forms first and possibly induces an eccentricity, but subsequent mass transfer 
from the white-dwarf progenitor is expected to circularize the orbit (Verbunt 
k Phinney 1995). Indeed PSR J0751+1807 (Lundgren, Zepka, k Cordes 1995; 
Nice, Splaver, k Stairs 2004) and PSR J1757-5322 (Edwards k Bailes 2001) 
have very small eccentricities (e < 10 -4 ). However, the discovery of the white 
dwarf companion to PSR B2303+46 (van Kerkwijk k Kulkarni 1999) led various 
groups to consider a different evolutionary path that could explain the observed 
eccentricity of e = 0.658 (Portgies Zwart k Yungelson 1999; Tauris k Sennels 
2000; Nelemans, Yungelson, k Portegies Zwart 2001; Brown et al. 2001; Davies, 
Ritter, k King 2002). 

In these modified scenarios, the original primary star, which is massive 
enough to form a white dwarf but not a neutron star, evolves first and transfers 
its mass to the secondary star. After the primary star loses its envelope it 
becomes a white dwarf, while the secondary continues its evolution but at a 
higher mass than its initial mass (high enough to form a neutron star) due to 
the mass transfer phase from the original primary. The primary eventually fills 
its Roche lobe and a common envelope phase ensues: the white dwarf and the 
helium core of the secondary spiral together as the common envelope of the two 
stars is ejected from the system. The outcome is a tight binary with a white 
dwarf and a helium star that is massive enough to explode in a supernova. This 
introduces an eccentricity to the orbit of the NS— WD binary that forms. Hence 
eccentric NS— WD binaries are formed because of two main facts relevant to 
the early evolution of the binary progenitor: (1) the initial binary components 
are both massive enough to form only a white dwarf, but still close to being 
massive enough to form a neutron star; (2) the initial mass transfer phase from 
the primary to the secondary increases the mass of the secondary enough that it 
can now form a neutron star eventually, but only after the initial primary ends 
its evolution as a white dwarf. 
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Table 1. Observational properties of eccentric NS-WD binaries. The columns 
indicate the pulsar name, spin period P, spin-down rate P, orbital period Pi,, the 
estimated mass of the WD companion m c , orbital eccentricity e, characteristic age 
t c , time to reach the death line t&. 
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3. Empirical NS— WD formation rate estimates 

In KKL, we introduced a statistical method to calculate a PDF of the rate 
estimates for Galactic close DNS systems. This method can be applied to any 
type of pulsar population of interest. 

Here, we consider eccentric NS— WD binaries. In order to calculate their 
Galactic formation rate TZ\>, we need to estimate: (1) the number N to t of Galactic 
pulsar populations with pulse and orbital characteristics similar to those in the 
observed sample of eccentric NS-WD (i.e. J1141-6545 or B2303+46); (2) the 
lifetime Tnf e of each observed system; (3) an upward correction factor for 
pulsar beaming. 

We calculate A?tot by modeling in detail the pulsar survey selection effects 
associated with the discovery of these systems for a number of parent pulsar 
population models described in KKL. The model assumptions for the pulsar 
luminosity function dominate the systematic uncertainties of our overall calcu- 
lation. 

For our reference model (model 6 in KKL), we obtain the most likely value 
of Ntot for PSRs J1141-6545 and B2303+46 to be 7V 1141 -370 and N 23 03 -240, 
respectively. 

The lifetime of the system is defined as ruf e = r c + r<j, where r c is the 
characteristic age. Td is the time which a pulsar will reach the "death line" (Ru- 
derman <fe Sutherland 1975). We calculated Td following (Chen & Ruderman 
1993; eq. (9) in their paper). The calculated PDF of the formation rate is dom- 
inated by PSR J1141-6545-like population due to both the shorter lifetime and 
number abundance of this population (Fig. 1). The estimated formation rate 
for our reference model is TZb= Q-8~^1'1~^1 3 q 7 Myr -1 at 68% and 95% confidence 
limits, respectively. The most likely values of TZ^ lie in the range 0.5 — 16 Myr -1 
for all the PSR population models we consider (see Table 2 for the rate estimates 
for a number of models with different pulsar luminosity functions). In the ab- 
sence of observational constraints on the geometry of both pulsars, we decided 
to not include any such upward correction for pulsar-beaming. Therefore, our 
estimated rates should be considered as lower limits. We can also compare these 
estimates with those for the Galactic DNS rate. The most likely values of DNS 
rates are found in the range 4 — 224 Myr -1 for the models considered with pulsar 
beaming correction, (see contribution by Kim et al. in these proceedings and 



4 V. Kalogera, C. Kim, D. R. Lorimer, M. Ihm <fc K. Belczynski 



5 



_ 112303 + 46 














66% - 


















95% 






\ 




V 99% 


" J 1 14 1 6545// 








\ 
\ 


\ 
\ 
\ 


\ 

i \ 


9xl0 4 


















6x10* 




68% 














"3x10* 




\s>f>% 
















■ 
















10 

1 


20 


:to 


i 






..i..., 1 



Galactic formation rate of ecc. NS— WDs (Myr 1 ) 

Figure 1. The PDF of Galactic formation rate estimates for eccen- 
tric NS— WD binaries (solid line) for our reference model. Dashed 
and dot-dashed lines represent the individual PDFs of the formation 
rates for sub-population of binaries similar to either PSR B2303+46 or 
J1141— 6545. No corrections for pulsar beaming have been applied. 



Kalogera et al. (2004)). Since the beaming correction factor for young pulsars is 
presumably larger than those for recycled pulsars, the true Galactic formation 
rate for eccentric NS— WD binaries are expected to be comparable to that of 
DNS systems we estimated. 

Population synthesis calculations have been widely used to study the forma- 
tion of various types of compact object binaries including eccentric NS— WD bi- 
naries. The range of details in the evolutionary calculations as well as the extent 
(if any) of the parameter studies vary significantly. In this section we summarize 
the main results from theoretical studies in literature and compare them to the 
empirical rates derived in the previous section. 

Portegies Zwart & Yungelson (1999) obtained a formation rate for eccentric 
NS— WD binaries comparable, but somewhat larger than that of DNS systems 
(7£b=44 and 34 Myr -1 for eccentric NS— WD and DNS systems respectively). 
Tauris & Sennels (2000) presented similar results (TZb^ 57 Myr -1 ) and con- 
cluded that the formation rate of eccentric NS— WD binaries is ~18 times higher 
than that of DNS systems. Davies, Ritter, & King (2002) estimated the Galac- 
tic formation rate of systems like J1141— 6545 or B2303+46. They obtained 
formation rates in the wide range of ~ 10~ 5 — 10~ 3 yr" 1 when considering both 
J1141— 6545— like and B2303+46— like systems with different evolutionary his- 
tories. Finally, Nelemans, Yungelson, & Portegies Zwart (2001) derived 240 
Myr -1 for the formation rate of all types of NS— WD binaries, which is compa- 
rable to other studies. We consider this rate to be an upper limit of eccentric 
NS— WD (shown as an open square with a downward arrow in Fig. 2). 

In addition to the above results from then literature, we have performed 
population synthesis calculations using the StarTrack population synthesis code 
(Belczynski, Kalogera, Sz Bulik 2002; Belczynski et al. 2004). We have explored 
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Table 2. The estimated Glactic formation rate and the most likely value of TVtot 
of eccentric NS-WD binaries for models with different pulsar luminosity functions. 
Model number is same with KKL. We show the most likely value of Rb at 68% and 
95% confidence limits. 
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a selection of models (more than any of the other studies) that promise to give 
us the widest variations of rate estimates. We derive a range ~ 30 — 100 Myr -1 
(using an absolute normalization of models based on empirical supernova rate 
estimates for our Galaxy; Cappellaro, Evans, & Turatto 1999). 

The different results from various population studies are mainly attributed 
to varying assumptions about the initial mass function, the initial mass-transfer 
process, the assumed star-formation or supernova rate used as a normalization 
factor, the initial binary fraction, and NS kicks. We also note that the lifetime 
of the system is a free parameter in the theoretical rate estimation, which is 
attributed to at least an order of magnitude of the uncertainty in the calculation. 



4. Comparison with rates from binary evolution 

In Fig. 2, we overlap the empirical formation rate estimates for eccentric NS— WD bi- 
naries (filled circles with error bars labed by KKL) with previous studied (open 
squares and a think solid line labeled by N01 (Nelemans, Yungelson, & Porte- 
gies Zwart 2001), TS00 (Tauris, & Sennels 2000), and PZY99 (Portgies Zwart & 
Yungelson 1999), and D02 (Davies, Ritter, Sz King 2002), respectively) as well as 
results from StarTrack (thin solid lines with filled triangles). It is encouraging 
that the overall rate estimates from different methods appear to be consistent 
with one another. If we consider the most likely values of the empirical rate 
estimates (filled circles), they are somewhat smaller than the theoretical values. 
However, at a 95% confidence interval, the empirical and theoretical estimates 
become comparable. We also emphasize that the upward correction for pulsar- 
beaming has not been applied, and therefore the empirical rate estimates should 
be considered as lower limits. 
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Figure 2. Comparison between the empirical and theoretical rate 
estimates. Error bars with filled triangles indicate results from 
Star Track, open squares and a solid line are adapted from the lit- 
erature, and filled circles with error bars are obtained in this work (see 
text for details). 



In conclusion, we find that our empirical estimates for the Galactic forma- 
tion rate of eccentric NS— WD binaries are overall consistent with the estimates 
derived based on binary population synthesis models. Despite the large uncer- 
tainties we consider this consistency as evidence that our current theoretical 
understanding for the formation of eccentric NS— WD is reasonable. However, 
the extent of the range covered by the empirical estimates and the population 
synthesis studies that attempt even a minimal parameter study (e.g., Davies, 
Ritter, & King (2002) and our results from StarTrack) indicates once again 
the necessity for careful parameter studies of rate calculations. It also indicates 
that the empirical rates could in principle be used to constrain binary evolution 
calculations. 
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